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ABSTRACT 



Context. We present the results from a study of the long-term optical spectral variations of BL Lacertae, using the long and well- 
sampled B and i?-band light curves of the Whole Earth Blazar Telescope (WEBT) collaboration, binned on time intervals of 1 day. 
Aims. To study the relation between the long-term spectral variations and the respective flux variations of the source. 
Methods. We use cross-correlation techniques to investigate whether there are any delays between the flux variations in different 
energy bands and between the flux and spectral variations. 

Results. The relation between spectral slope and flux (the spectrum gets bluer as the source flux increases) is well described by a 
power-law model, although there is significant scatter around the best-fitting model line. To some extent, this is due to the spectral 
evolution of the source (along well-defined loop-like structures) during low-amplitude events, which are superimposed on the major 
optical flares, and evolve on time scales of a few days. The B and R-hand variations are well correlated, with no significant, measurable 
delays larger than a few days. On the other hand, we find that the spectral variations lead those in the flux light curves by ~ 4 days. 
Finally, during at least the largest amplitude flares, the S-band variations appear to evolve faster than those in the R band. 
Conclusions. We confirm the "bluer-when-brighter" mild chromatism of the long-term variations, and we show that it can be explained 
if the flux increases/decreases faster in the B than in the R band. We also report the discovery of the lag between spectral and flux 
changes. These two features can be explained in terms of Doppler factor variations due to changes in the viewing angle of a curved 
and inhomogeneous emitting jet. 
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1. Introduction 

BL Lacertae is the prototype of a class of active galactic nuclei 
(AGNs) known as "BL Lac objects", or simply "BL Lacs". The 
members of this class show non-thermal continuum energy spec- 
tra, fast and large-amplitude variations from radio up to y-rays, a 
high degree of linear polarization, and radio jets with individual 
components often exhibiting apparent superluminal motion. 

, The "Whole Earth Blazar T elescope" (WEBT; 

' | http : //www ■ to ■ astro . it/blazars/webt/) is a large in- 
ternational collaboration among optical and radio astronomers, 
established in 1997. It organizes monitoring campaigns on 
selected blazars in order to obtain continuous, high-temporal- 
density light curves in the optical and radio bands (usually in 
conjunction with o bservations at other wavelengths like X and 
y-rays; se e e.g. Bottcher et al. 2005; Raiteri e t alj|2005l 120061: 
lOstorero et al. 2006 ; Villata et al, 2006, 2007). 

BL Lacertae has been the target of four WEBT monitor- 
ing campaigns in the past. The first two were rather short. 
They were organized in 19 99, simultaneousl y with ASCA and 
BeppoSAX observations ( Ravasio et al.l 120021) . The remaining 
two were long-term ca mpaigns cair ied out in the p eriods May 
2000 - January 2001 (iVillata et alJ j2002; Bottcher et al.ll2003t) 



and May 2001 - February 2002 (iVillata et al.ii2004bllal) . Both 
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campaigns resulted in long, well-sampled, and high-precision 
light curves. 

Using data fro m the first long-term WEBT campaign, 
IVillata et all (|2002| ) found that the B and /J-band light curves 
were well correlated with no measurable time delay. They also 
found that the flux variations are associated with significant 
spectral variations. They interpreted them in the context of a two 
variability mechanism model: the first mechanism is essentially 
achromatic and it is responsible for the long-term variations (i.e. 
variations which operate on time scales of a few days), while the 
second one causes fast (i.e. shorter than a day) flares, superposed 
on the long-term variations, and introduces spectral changes in 
the sense that the spectrum becomes harder (bluer) as the source 
brightens. 

Similar results were obtained bv IVillata et al.l (l2004bl) . who 
studied not only the light curves from the 2001-2002 WEBT 
campaign, but also composed light curves from 1994 to 2002, in 
all optical bands, with data taken from the literature as well as 
by members of the WEBT. The authors distinguished a "mildly- 
chromatic" (in the sense that the B - R versus R plot has a slope 
~ 0. 1), long-term variability component, which operates on time 
scales longer than a few days, from a "strongly-chromatic" (in 
which case the variations trace a 0.4 slope in the B - R versus R 
plot), short-term variability mechanism, which operates on intra- 
day time scales. 

IVillata et al.l ( l2004bh based their results mainly on the study 
of the best-sampled parts of the light curves, which cover the 
period between 1997 and 2002. In this work, we focus our at- 
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tention on the B and R-hand light curves of BL Lacertae in the 
same period. Our main aim is to better understand the long-term 
chromatic behaviour of the source, and to investigate possible 
reasons for the observed spectral variations. 

In Sect. 2 we present the light curves we use in this work 
and we quantify the "spectral slope versus flux" relation on time 
scales longer than a day. Sect. 3 deals with the cross-correlation 
analysis between the B and /?-band light curves, and between the 
spectral slope and flux variations. In Sect. 4 we present the re- 
sults from the study of two long, well sampled, large-amplitude 
flares in the light curves, while our final conclusions are dis- 
cussed in Sect. 5. 

2. The optical light curves 

IViflata et al] ( l2004bl) presented UBVRI Ught curves of 
BL Lacertae from 1994 up to 2002, using data from the four 
WEBT campaigns and the literature. The host galaxy constant 
flux contribution was subtracted from the observed light curves 
in all bands as explained in Section 5 of Villata et al., 2002. We 
chose to use the Villata et al. (2004b) B and /?-band light curves 
from the period between 1997 and February 2002, as we be- 
lieve that these are the "best" light curves (in terms of length, 
number of data points and sampling frequency of observations) 
among the currently available optical datasets of this object. The 
data that we use in this work are now publicly available in the 
WEBT archive (http://www.to.astro.it/blazars/webt/) and can be 
requested by sending a message to the WEBT president. 

In order to minimize any systematic ef fects that cou l d influ- 
ence our results we examined carefully the lViUata etal] (l2004bh 
data and selected our final B and 7?-band datasets as follows: 
a) we kept data with errors not greater than 0.04 and 0.03 mag 
in the B and /?-band light curves, respectively; b) we used only 
those data points which have at least one counterpart in the other 
band within 20 min, either from the same telescope or from the 
same literature paper; c) we also removed three "unreliable" B 
data points. 

The two upper panels in Fig. [T]show the B and /?-band light 
curves (in flux density units) that we use in this work. Since we 
are interested in studying the spectral variations on time scales 
longer than a day, the light curves are binned in 1 -day long in- 
tervals, starting from JD = 2450644. If more than one point is 
present in the interval, we take their average flux density and 
time of observation. We consider the standard deviation, i.e. Icr, 
the average spread of the points around their mean, as the "error" 
on the average flux density in the respective bin. In the case there 
is just one point in a bin we accept its error as it is. Typically, 
the error of the points, estimated in this way, in the B and R 
band 1 day binned light curves is of the order of ~ 0.1 - 1 mJy. 
Therefore, although we do plot the error bars of the points in 
Fig. [T] most of them are smaller than the size of the symbols we 
use. 

The source is brighter in the R band. The maximum flux den- 
sity recorded in that band is almost twice the respective maxi- 
mum flux density in the B-band light curve. Apart from this dif- 
ference, the two light curves look very similar The same varia- 
tions appear, simultaneously, in both of them. However, the am- 
plitude of the observed variations is larger in the B-band light 
curve. Indeed, the "fractional variability amplitude", /,n,s (de- 
fined as th e ratio of the standard deviation over the light curve 
mean, see iPapadakis et al.l l2003l) . of the B and R light curves 
is 65.3 + 0.2% and 57.1 + 0.2%, respectively. The errors of 
/nns account only for the measurement errors in the light curve 
points, and have been estimated according to the prescription of 
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Fig. 1. First and second panel from top: The 1997-2002 B and R 
light curves of BL Lacertae binned in intervals of size equal to 1 
day. Third panel from top: The Fb/Fr flux density ratio, plotted 
as a function of time. Bottom panel: Fb/Fr plotted as a function 
of (Fb + Fr)/2. The solid line shows the best fitting power- law 
model to the data. 



IVaughan et alj|2003 l. Given the very good correlation between 
the B and R band light curves, one could assume that the ~ 1% 
difference in frms,B and /rms.R B maybe due to differences in 
the experimental noise level of the two light curves. For ex- 
ample, Poisson noise will introduce some scatter in the /rms 
values we measure. However, the errors we provide (following 
Vaughan et al. 2003) suggest that this is not the case. The differ- 
ence between the B and R band variability amplitudes is most 
probably an intrinsic property if the variability mechanism that 
opearates in the source. 

2.1. Spectral variations 

The third panel from the top in Fig.[T]shows a plot of the Fb/Fr 
flux density ratio as a function of time. This ratio can be consid- 
ered as representative of the slope of the optical spectrum of the 
source. Errors on this ratio have been estimated using th e usual 
"propagation of errors" recipe (see e.g. lBevingtonll 19691) . 

Since the B and /?-band light curves do not have the same 
fractional variability amplitude, we expect the spectral slope 
(and hence the ratio Fb/Fr) to change with time. Indeed, 
Fig. [T] shows clearly that this ratio is highly variable. The frac- 
tional amplitude, /rms.R, of the Fb/Fr variations is 10.3 + 0.2% 
(where the error accounts only for the measurement error of the 
Fb/Fr points, as explained above). We conclude that, on time 
scales longer than 1 day, the optical flux density variations of 
BL Lacertae are associated with moderate, but significant, spec- 
tral slope variations as well. 

In the bottom panel of Fig. [T] we plot the Fb/Fr ratio as 
a function of (Fb + Fr)/2, i.e. the mean of the B and 7?-band 



Papadakis et al.: Optical spe 



flux densities, which can be considered as representative of the 
flux density in an intermediate band. We have chosen to use 
this quantity in the study of the "spectral slope (Fb/Fr) ver- 
sus brightness" relation, as the average of the source signal in 
the two bands should minimize the possibility of introducing ar- 
tificial correlations in this relation. Such correlations could be 
introduced if we were using either the B or the /?-band flux den- 
sities, individually, because of the interdependency of the Fb/Fr 
values nd these measurements. 

The "spectral slope vs brightness" plot shows clearly that 
the flux and spectral variations are well correlated: as the source 
flux increases (decreases), the spectrum becomes harder (softer). 
The solid line in the bottom panel of Fig. [T] shows the best 
fitting power-law (PL) model to the data: Fb/Fr oc [(Fb + 
Fr)/!]^-^'^'^'^'^^ . Although it does not provide a statistically ac- 
cepted fit (we get ax^ value of 6313.5 for 308 degrees of free- 
dom when we consider the measurement uncertainty of each 
point in the plot), it nevertheless describes rather well the overall 
trend of the Fb/Fr increase with increasing flux. 

The statistically poor quality of the fit is due to the fact 
that there exists a significant scatter in the Fb/Fr values for a 
given flux density. However, it is of small amplitude (at each 
flux level, the maximum deviation around (/^B/f'R)niodei is less 
than ~ +15% of {FB/FR)mode{) and does not show any systematic 
deviations from the best fitting PL model. This scatter implies 
some residual, weak spectral variability which is not correlated 
with the flux density level of the source. 

Interestingly, the Fb/Fr oc [(Fb + Fr)/2]^-^'^ relation would 
correspond to a line of slope ~ 0.14 in a "B - /? vs magnitude of 
(Fb + Fr)/2" plot. This is close to the "s ofter" slope of 0. 1 in the 
"B -Rws R" plot of lVillataetaLl (l2004 b). Hence the flux-density 
variations that we observe in the 1-day binned light curves of 
Fig. [1] corres pond to the "mildly-chromatic", long-term varia- 
tions of (.ViUata et al.ll2004h: see also .Hu et al..,2006i) . 



3. Cross-correlation analysis 

3.1. Correlation between the B and R -band light curves 

A possible explanation for the observed spectral variations is the 
existence of a delay between the B and 7?-band light curves: 
if the flux increases or decreases with the same rate in both 
light curves, but the B-band variations lead those in the R 
band, then the spectrum is "bluer"/"redder" when the flux is ris- 
ing/decaying, respectively. 

In order to investigate this possi bihty we used the discret e 
correlation function (DCF) method of lEdelson & Krolikl(ll988l) . 
Our results, with a DCF bin size of 4 days, are plotted in the 
top panel of Fig. |2](open circles). We get similar results when 
we use a wide range of DCF bin sizes, from 1 up to 10 days. A 
positive lag of the DCF peak in this plot would imply that the 
Z?-band variations lead those in the R band. In the same panel we 
also plot the B-band auto-correlation function (ACF; solid line). 
The ACF of the Fr and (Fb + Fr)/2 light curves are practically 
identical to that of the Fb light curve. 

The Fb vs Fr DCF is similar to the B-band ACF and shows 
a strong (DCFn,ax ~ l)and narrow peak at zero time lag. Strong 
and narrow DCF maxima at zero lags can be introduced when 
systematic errors affect, in the same way, the data points in both 
light curves. However, in our case, the points in the light curves 
shown in Fig. [1] correspond to observations from many different 
telescopes, so that the presence of a global systematic effect is 
rather unlikely. 
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Fig. 2. Top panel: The B-band ACF (solid line) and the Fb vs 
Fr DCF (open circles). Bottom panel: The (Fb + Fr)/2 and Fb 
vs Fb/Fr DCFs (open circles and solid line, respectively). The 
positive-lag Fb vs Fb/Fr DCF part is reflected about the zero- 
lag axis (dashed line) for comparison reasons (see text for de- 
tails). 



If we omit the DCF point at zero lag, then the highest DCF 
value is ~ 0.7. If we now estimate the centroid lag, Tcent, using 
all DCF points with values in excess of 0.7 DCF^ax (assuming 
DCFmax ~ 0.7), we find that Tcent - 2 days. This result suggests 
that there may exist after all a small delay between the variations 
in the B and /?-band light curves. The positive centroid lag is due 
to the fact that at negative lags (down to ~ -60 days) the Fb vs 
Fr DCF is smaller than the Fb ACF while the opposite trend, i.e. 
Fb vs Fr DCF being slightly larger than the Fb ACF, is observed 
at the positive lags up to 60-70 days. Admittedly, this is a low 
amplitude effect, and longer light curves are necessary in order 
to test its validity. 

3.2. Correlation between flux and spectral variations 

We also investigated the cross-correlation between the flux- 
density light curves and the Fb/Fr curve. In the bottom panel 
of Fig. |2] we plot the (Fb + Fr)/2 vs Fb/Fr and the Fb vs 
Fb/Fr DCFs (open circles and solid line, respectively; the Fr 
vs Fb/Fr DCF is very similar to the ones plotted in this panel). 
Positive lags mean that variations in the flux-density light curve 
lead those in the flux ratio curve. 

The two DCFs look very similar. They show a broad hump 
around zero time lag, with DCF^ax ~ 0.7 in both cases. This 
result suggests that the flux and spectral variations are well cor- 
related. This is not surprising, given the Fb/Fr vs (Fr + Fr)/2 
plot shown in the bottom panel of Fig. [1] The fact that the max- 
imum DCF value is not larger than 0.7 can be explained by the 
significant, albeit of low amplitude, scatter of the points around 
the best fitting PL model in this plot. 
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Fig. 3. The cumulative distribution function of Tcent in the case 
of the (Fb + Fii)/2 vs Fb/Fr DCF, as estimated using the resuhs 
of the numerical experiment that we describe in Sect. 3.2. 



The most interesting aspect of both DCFs is the strong asym- 
metry towards negative lags. In order to highlight this effect, in 
the bottom panel of Fig. |2] we also plot the positive-lag Fb vs 
Fb/Fr DCF curve reflected about the zero-lag axis (dashed line). 
It is now evident that the amplitude of the DCF at lags smaller 
than ~ - 10 days is larger than that at the respective positive lags. 
This implies that there are delays between the flux and spectral 
variations, with the latter leading the former 

In order to quantify the average delay between the spectral 
and flux variations, we estimated the centroid of the (Fb + Fr)I2 
vs Fb/Fr DCF using the same method as above. We found 
that Tcent - -4 days. In order to estimate the uncertainty in 
this measurement, w e employed the bootstrap techniques of 
iPeterson etalldigQSl) and created 10000 pairs of simulated light 
curves. In Fig. [3] we plot the sample cumulative distribution of 
the centroid value as was estimated using the Tcent values that we 
computed for the 10000 synthetic DCFs of our numerical exper- 
iment. 

The average centroid value is equal to -4 days, and 68.3% 
of the synthetic DCFs yield centroid values between -6.4 and 
-2 days. Furthermore, only 2.4% of all DCFs resulted in Tcent 
having a positive value. We conclude that the delay of ~ 4 days 
that we detect between the Fb/Fr and the flux variations is sig- 
nificant. 

We should point out that the determination of the DCF 
centroid is based on a rather subjective method and depends 
on the number of the DCF points involved in the calculation. 
Furthermore, the DCFs shown in the bottom panel of Fig.|2]show 
broad humps rather than well defined, narrow peaks. The broad- 
ness of the DCF maxima may imply the presence of more than 
just a single delay between the spectral slope and flux variations. 

Having these remarks in mind, we conclude that, although 
the flux ratio Fb/Fr follows rather well the increase/decrease of 
the source flux (see bottom panel in Fig.[T]), the cross-correlation 
analysis reveals a subtle detail: Fb/Fr decreases or increases 
before the respective flux changes. This result is significant at 
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Fig. 4. Top panel: The last part of the B and /?-band light curves 
of the 1997 flare. The data are normalized to the flux density of 
the rightmost point. Middle panel: Fb/Fr plotted as a function 
of time, using the data shown in the upper panel. The vertical 
line indicates the time of maximum flux. Some points are plotted 
with asterisks and filled squares for reasons explained in the text. 
Bottom panel: Fb/Fr plotted as a function of (Fb + Fr)/2. The 
solid line shows the best-fitting PL model to the data. 



the 97.6% level. Our best estimate of the average delay between 
the flux and spectral slope variations is -4;^2 4 days. 

In order to further investigate this interesting issue, as well 
as the origin of the spectral variations, in the following section 
we focus our attention on two best sampled flares in the optical 
light curves shown in Fig.[T] 

4. Spectral analysis of individual flares 

4.1. The 1997 flare 

The 1997 and 2001 large-amplitude optical flares of BL Lacertae 
are well sampled in the light curves that we use in this work (the 
respective data points in Fig. [T] are those in the time intervals 
5650-5750 and 7050-7150). 

In the top panel of Fig. |4] we plot the B and /?-band light 
curves corresponding to the last part of the 1997 flare (filled and 
open circles, respectively). The data sample the last few days of 
the high-brightness phase and the subsequent flux decline phase 
which lasted for ~ 50 days. We do not plot the data before JD 
- 2450680 because they are erratic and do not show any other 
well-defined flux rising or decaying parts. The light curves are 
normalized (F'^, F'^^) to their rightmost point (which, as a result, 
is set equal to one, for both of them), as this is their lowest flux 
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point. In this way we can better compare variability amplitudes 
between the two light curves. 

At the flare peak, the B and 7?-band flux densities are 21.1 + 
0.6 and 15.6 + 0.4 times higher than the respective lowest lev- 
els. The diff'erence between these two values (which represent 
the maximum variability amplitude of the light curves shown in 
Fig. nil is 5.5 + 0.7. However, as the flux decreases, this differ- 
ence decreases as well. For example, ~ 35 days after the flare 
peak, the B and /?-band flux densities are only 1.45 + 0.04 and 
1.41 + 0.02 times larger than the respective lowest points. From 
then on, the B and /?-band variability amplitudes are quite sim- 
ilar, indicating that the flux in the B band decreased /oifer than 
the flux in the R band. 

In the middle panel of Fig.|4]we show the Fb/Fr ratio, which 
is clearly variable. This result implies that the optical spectrum 
of the source changes, becoming "softer" ("redder") as the flux 
decreases. These spectral variations are due to the akeady no- 
ticed faster decrease of the B-band flux. 

Interestingly, the data show that Fb/Fr was decreasing even 
before the flux reached its maximum level. The vertical, dashed 
line in the top panels of Fig. |4] indicates the time of the flare 
peak. Although the errors of the two previous FbIFr points are 
large, the data indicate that the flux ratio started decreasing in 
advance of the source flux. This behaviour is consistent with the 
"{Fb + Fk)/2 vs Fb/Fr" DCF results that we reported in the 
previous section. 

The bottom panel of Fig. |4] shows the "Fb/Fr vs (Fb + 
Fr)I2" plot using the data that are plotted in the top panels of 
the same figure. The solid line shows the best-fitting PL model 
to the data. Its slope is ~ 0.13, almost identical to that of the 
best-fitting PL model to the data shown in the bottom panel of 
Fig.[T] We conclude that a constant difference between the vari- 
ation rates in the two bands can produce spectral variations sim- 
ilar to the long-term, "mildly-chromatic" spectral variations that 
we observe in BL Lacertae between 1997 and 2002. 

In the middle and bottom panels of Fig. |4] we have used dif- 
ferent symbols for some points which appear to deviate signif- 
icantly from the generally "smooth" trend of Fb/Fr decreasing 
with time. Points plotted with asterisks and filled squares corre- 
spond to the flux "humps" that we observe in the time intervals 
5707-57 12 and 57 12-5722, respectively. Although the flux vari- 
abiUty amplitude of these "mini-events" is small, the associated 
spectral variations introduce significant scatter around the best- 
fitting PL model in the ''Fb/Fr vs (Fb + Fr)/2" plot. 

We may thus speculate that some of the scatter around the 
best-fitting PL model to the data shown in the bottom panel of 
Fig. [1] is introduced by short, low-amplitude events which are 
accompanied by spectral variations that do not follow the general 
trend of the larger-amplitude flares. 

We note that there exists an interesting correlation between 
the spectral variations of these low-amplitude events and the 
respective flux changes. We have connected with arrows the 
points plotted with asterisks and filled squares in the "Fb/Fr vs 
{Fb + Fr)!!" plot of Fig.|4]to indicate the time evolution of these 
"mini-events" in this plane. We can see that these points define 
loop-like structures which evolve in the anti-clockwise direction. 

4.2. The 2001 flare 

During the BL Lac 2001 WEBT campaign, two large-amplitude 
flares, the first in mid 2001 and the second in early 2002, 
were detected ( Villata et al. 2004b). The 2001 flare was double- 
peaked. The second peak is well sampled in the WEBT light 
curves that we use in this work. The respective data are plotted 
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Fig. 5. Same as Fig. H] but for the 2001 flare. The solid and 
dashed lines in the top panel show the best-fitting exponential 
functions to the B and 7?-band data, respectively. The solid line 
in the middle panel shows the resulting Fb/Fr ratio. Filled and 
open squares in the bottom panel show the data during the flux 
rising and decaying phases, respectively. Some points in the mid- 
dle and bottom panels are plotted with asterisks for reasons ex- 
plained in the text. 



in the top panel of Fig. |5] Just like in Fig. |4] we have normalized 
both light curves to the flux density of their rightmost point. In 
the middle panel we show Fb/Fr as a function of time, and in 
the bottom panel the same quantity is plotted as a function of 
(Fb + Fr)/2. The advantage of the 2001 light curves is that they 
sample well both the flux rising and flux decaying parts of the 
flare. 

We fitted the rising and decaying parts of the light curves 
with an exponentially increasing, F(t) oc exp[(f - f())/ri], and 
decreasing, F(t) oc exp[-(f - fo)/T"d], function, respectively (fo 
represents the time of the flare peak). The best-fitting models are 
also plotted (normalized as the light curves) in the top panel of 
Fig. |5] (solid and dashed line for the B and /?-band light curves, 
respectively). 

The model describes rather well the overall flux evolution, 
during both the rising and decreasing phases of the flare. Our 
best-fitting t values are: = 16,rr.R - 18,Td.B = 13, and 
Td,R = 15 days (with an error of +0.8 days in all cases). The 
differences among these best-fitting r values suggest that: a) the 
flare is not symmetric, as the flux decrease is faster than the flux 
increase, i.e. is smaller than in both bands (by the same 
amount of ~ 3 days), and b) the flare evolves faster in the B than 
in the R band, as both the rise and decay tb values are smaller 
than the respective tr values by about 2 days. 
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The second result can explain the observed spectral varia- 
tions. The solid line in the middle panel of Fig. |5] shows the 
Fb/Fr curve derived from the best-fitting B and /?-band expo- 
nential functions. Obviously, the Fb/Fr model curve fits the data 
well. At the beginning of the flare, the spectrum hardens as the 
flux increases, because the Z?-band flux increases faster than the 
flux in the R band (ti-^r > Ti-b)- After the flare peak, the spec- 
trum softens, because the B-band flux decreases faster than the 
/?-band flux. 

The vertical dashed line in the upper two panels of Fig.|5]in- 
dicates the time of the flare peak. One can see that the change 
of the spectral evolution from spectral hardening (i.e. Fb/Fr 
increasing) to spectral softening (Fb/Fr decreasing) happens 
a couple of days before the flux reaches its maximum. This 
is consistent with the DCF results that we report in Sect. 3.2. 
Furthermore, the data points in the "spectral slope vs time" plot 
indicate that the change in the spectral evolution is rather abrupt. 
This is also clearly indicated by the significant "break" that ap- 
pears in the Fb/Fr model curve near the flux peak. 

In the bottom panel of Fig. |5] we show the "Fb/Fr vs 
(Fb + Fr)/2" plot. Filled and open squares correspond to the 
data during the rising and decaying phases, respectively. A PL 
model of slope ~ 0.13 (not plotted for clarity reasons), equal to 
that of the 1997 flai-e, describes wefl the "Fb/Fr vs {Fb + Fr)/T 
relation. 

The solid (blue) lines in the same panel show the "spectral- 
slope vs brightness" relation when we consider the best-fitting 
exponential rising and decaying functions, which are plotted in 
the upper panel of Fig. |5] In agreement with the data points 
(open squares lay systematically below the points plotted with 
filled squares), and because of the abrupt change of the spectral 
slope near the flare maximum, the "spectral-slope vs brightness" 
model curve during the flux decay phase lies below the respec- 
tive flux rise model curve. 

Just like during the 1997 flare, low-amplitude flux variations 
introduce significant scatter in the "Fb/Fr vs (Fb + Fr)/2" plot. 
The points plotted with asterisks in the middle and bottom panels 
of Fig.|5]highlightthe spectral variability behaviour of the source 
during such "events". 

For example, during the time interval ~ 7093-7102 the flux 
decreased and then increased until the flare reached its peak. The 
respective points in the "spectral-slope vs flux" plot follow a well 
defined loop, which evolves clockwise (loop "I"), as the arrows 
clearly show in the bottom panel of Fig. |5] 

We observe loop-like patterns during the flux decaying part 
of the flare as well. The clearest example is the spectral evolu- 
tion during the end of the flare, i.e. after ~ 7117. can see that 
these points correspond smooth flux increase, which becomes 
more In the Fb/Fr vs (Fb + Fr)/2 plot, the respective data, plot- 
ted with asterisks, follow a clockwise evolving loop-like pattern 
(loop "11"). Note that Fb/Fr has already started increasing a few 
days before the flux increase. 



5. Discussion and conclusions 

BL Lacertae is one of the most frequently observed blazars, at all 
frequencies. It has been observed for more than a century in the 
optical band. The observations in the last years have been very 
intense, due to the recent WEBT campaigns to a large extent. 
The light curves that have resulted from these observations are 
densely sampled over a long time period (~ 8 years). Although 
optical observations of the source will be carried out in the future 
as well, we believe that the current light curves are of sufficient 



quality (in terms of length and sampling) to allow the investiga- 
tion of its"typicar' variability characteristics. 

Villata et al. (2004b) already showed that the variations in 
the optical light curves of BL Lacertae can be interpreted in 
terms of a "strongly" and a "mildly-chromatic" component, 
which operate on short (i.e. < 1 day) and longer time scales, re- 
spectively. They also showed that the short and longer time scale 
components determine a "bluer-when-brighter" slope of ~ 0.4 
and ~ 0.1, respectively, in a"B - R vs R" plot. 

In t his work, we continued the investigation of I Villata et al.l 
(I2004bh . and we studied in more detail the long-term, mildly- 
chromatic spectral variations that these authors had already no- 
ticed. Our main results are as follows: 

a) There is an indication that the /?-band variations on time 
scales longer than a day are slightly delayed with respect to the 
variations in the B band. However, this is a low amplitude effect 
and longer light curves are needed to confirm it. In any case, 
even if real, this effect cannot explain the significant, long-term, 
optical spectral variations that we observe. 

b) A constant difference between the rates with which the 
flux rises/decays in the B and /?-band light curves can explain 
the spectral variations that we observe in BL Lacertae. 

For example, we find that the 2001 flare data are consis- 
tent with exponentially rising/decaying functions. When the flux 
rises, Fb/Fr oc e\p[(t-tQ)iT,ji-Tr,B)/iTr,RTr,B)]- Using our best- 
fitting and Tr,B valucs, we find that Fb/Fr oc F^ '^. In the 
"B - R vs R" plane this relation corr esponds to a l ine of slope 
~ 0.12, consistent with the findings of IViUata etall (l2004bl) . We 
get a similar result when we consider the ratio Fb/Fr during the 
flux decay phase. In this case, Fb/Fr oc F^ '^. 

Obviously, the rise and decay time scales, as well as the max- 
imum variability amplitude, are not identical in all the "flare- 
like" events in BL Lacertae. In Fig.[T]one can detect many low- 
amplitude events, which last much less than the large-amplitude 
2001 flare. What our results suggest is that the characteristic 
time scales of the flares cannot change arbitrarily from one 
to the other. As an example, if we assume exponentially ris- 
ing/decaying flares, their time scales T,/d,R and Tr/d,B should 
change in such a way so that the ratio (rr/d,R - T^/d,B)/Tr/d,B al- 
ways remains roughly equal to ~ 0.1. In this way, the relation 
between the spectral slope and flux variations can be explained. 

c) We find that the flux variations follow the spectral vari- 
ations by an average delay of ~ 4 days. This is the first time 
that such an effect has been observed in BL Lacertae. A similar 
behaviour, in the optical band, h as also been observed recently 
in 3C 66A jBottcher et al.ll2005l) . This effect can be clearly seen 
in the well-sampled 1997 and 2001 flares. The DCF results we 
report in Sect. 3.2 suggest that this delay operates most of the 
times in BL Lacertae. 

If the "spectral-slope vs brightness" plot can be explained by 
a constant difference between the variability time scales in the 
various bands (as we showed above), then the spectral hysteresis 
effect can be explained if one assumes that the variability time 
scale does not remain constant during the flare flux rise phase. 
The spectrum can start softening before the source reaches its 
maximum flux if the variability time scale, which determines the 
flux evolution, starts increasing before the peak at a rate which 
increases with frequency, or vice versa. 

d) Although a power-law model describes well the "bluer 
when brighter" trend in the "spectral-slope vs brightness" plot, 
it cannot provide a statistically acceptable fit to it. There exists 
significant scatter around the best-fitting model line. 

To some extent, this scatter can be explained as the "rem- 
nant" of the fast, strongly-chromatic flares of IVillata et al.l 
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( l2004bl) in the 1-day binned light curves we use. It is quite possi- 
ble that in many cases the daily average points do not correspond 
to the "mean" of the fast flares but rather to a peak or a dip of 
the fast flares. 

At the same time though, as our analysis in Sect. 4 shows, 
some of the scatter is also due to the fact that, occasionally, the 
spectral evolution of the source defines loop-like structures in 
the "spectral-slope vs brightness" plot. In fact, if the ratio Fb/Fr 
decreases/increases a few days before Fb and Fr (see result "c" 
above), then most flares that last for more than 5-10 days should 
show clockwise evolving loop structures in the FbIFr vs flux 
plot. This is what we observe during the 2001 flare. 

Anti-clockwise evolving loops (like the ones observed in 
the flux decaying phase of the 1997 flare) imply that the spec- 
tral slope variations follow those in the flux light curves. The 
spectral-slope vs flux DCF results suggest that these must be 
rare events. In any case it will be interesting to investigate this 
issue in the future, using the even longer light curves that will be 
available then. 

In conclusion, our two main findings are: i) the confirma- 
tion of the bluer-when-brighter mild chromatism of the BL Lac 
long-term optical variations, which is not caused by any delays 
between the B and /?-band light curves, but can be explained if 
the flux rise/decay rates in the two bands are different, by a fixed 
amount, and ii) the discovery that the spectral-slope variations 
ten d to lead the flux vari ations, by a few days. 

IVillata et al.l (l2004bh explained the bluer-when-brighter be- 
haviour in terms of a Doppler factor, 6, variation on a convex 
spectrum. Indeed, an i ncrease in 6 imphes bo th an increasing 
flux {Fy ~ 6^; see e.g. IViUata & Raiterill 19991) and a change in 
the observed frequency (v ~ 6). The former effect causes the 
flux to increase at all frequencies, thus the appearance of a flare, 
while the intrinsic spectrum of the source is also "boosted" to- 
wards higher frequencies. If this spectrum is convex, we will ob- 
serve a harder and harder part of it in our optical band, during the 
rising part of the flare, and vice versa during the dimming phase, 
i.e. the bluer-when-brighter trend. In other words, a change of 
5 on a convex spectrum makes flux variations in B and R bands 
happen simultaneously, but the Z?-band flux rises/decays faster 
than the /?-band one. 

Moreover, if the 6 variation is due to a change of the jet 
viewing angle (<5 - [yil - /jco sff)]"') and the jet is inhomo- 
geneous and curved (see e.g. Villata et"al]|2007h . we will also 
have a tendency to a progressive reddening of the spectrum dur- 
ing the flare. Indeed, since the jet is inhomogeneous, it is com- 
posed of "slices" emitting different frequency ranges, which are 
redder and redder going outwards, along the jet curvature. Thus, 
as the jet approaches the line of sight, the 6 variation will af- 
fect first the bluer and then the redder parts, implying a "redden- 
ing" of the spectrum. However, this reddening is opposed by the 
"bluer-when-brighter" trend during the flare rising phase. The 
combination of the two effects can result in a softening of the 
spectrum before the flare peak, as observed. In other words, the 
concomitance of the "bluer-when-brighter" trend with the pro- 
gressive reddening, both of geometrical-kinematic origin, can (at 
least qualitatively) account for the changes in the spectral slope 
we analyzed in this paper. 

A further observing as well as modeling effort is needed 
to clarify and quantify this interpretation. On the observational 
side, apart from BL Lacertae, the WEBT archive akeady stores 
long and well sampled light curves for three more sources, and 
more data will appear in the near future. The analayis of these 
light curves in a way similar to what we presented in this work, 
and the comparison of the results with those presented here, will 



hopefully help us understand better the variability mechanism in 
BL Lacs. 
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